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In this work we study the dynamics of the axisymmetric Bianchi IX cosmological model with a term of 
quantum potential added. As it is well known this class of Bianchi IX models are homogeneous and anisotropic 
with two scale factors, A(t) and B(t), derived from the solution of Einstein's equation for General Relativity. 
The model we use in this work has a cosmological constant and the matter content is dust. To this model we 
add a quantum-inspired potential that is intended to represent short-range effects due to the general relativistic 
behavior of matter in small scales and play the role of a repulsive force near the singularity. We find that this 
potential restricts the dynamics of the model to positive values of A(t) and B(t) and alters some qualitative and 
quantitative characteristics of the dynamics studied previously by several authors. We make a complete analysis 
of the phase space of the model finding critical points, periodic orbits, stable/unstable manifolds using numerical 
techniques such as Poincare section, numerical continuation of orbits and numerical globalization of invariant 
manifolds. We compare the classical and the quantum models. Our main result is the existence of homoclinic 
crossings of the stable and unstable manifolds in the physically meaningful region of the phase space (where 
both A(t) and B(t) are positive), indicating chaotic escape to inflation and bouncing near the singularity. 
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I. INTRODUCTION 

Belinskii, Khalatnikov and Lifishitz yj started the 
question of chaotic behaviour of general Bianchi IX 
models in Relativistic Cosmology. The interest in 
the chaoticity (or not) of Bianchi IX models has 
been mainly focused on the Mixmaster case (vacuum 
Bianchi IX models with three scale factors ]2|])- The 
question of the generic behaviour (chaotic or not) of the 
Mixmaster dynamics remained unsettled mainly due to 
the absence of an invariant (or topological) character- 
ization of chaos in the model (standard chaotic indi- 
cators as Liapunov exponents being coordinate depen- 
dent and therefore questionable 13l|4[]). For discussions 
of the issue of chaotic dynamics on these models we re- 
fer to the works of lH-ll. 

In the early literature Hawking and Page lUOIl have 
shown that Bianchi IX models with a scalar field 
present a fractal set of uncountably infinite oscilla- 
tory orbits, which could be considered a signature of 
chaos in the dynamics. Cornish and Levin [11] pro- 
posed to quantify chaos in the Mixmaster universe 
by calculating the dimensions of fractal basin bound- 
aries in initial-conditions sets for the full dynamics, 
these boundaries being defined by a code association 
with one of the three outcomes on which one of the 
three axes is collapsing most quickly, as established 
numerically. Joras and Stuchi I112ll examined chaos in 
FRW models with a coupled scalar field by extending 
the analysis to the complex plane and found a family 
of non-collapsing structures. In particular homoclinic 
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chaos in axisymmetric Bianchi IX universes with mat- 
ter and cosmological constant has been treated in IU3I1 
using coordinate independent topological structure of 
the dynamics to characterize chaotic behavior in the 
Hamiltonian system and its physical implications. This 
paper follows along these lines. 

The phase space of the classical model is noncom- 
pact and the presence of the cosmological constant de- 
termines two crucial facts in phase space: first, the 
existence of a critical point of the saddle-center type; 
second, two critical points at infinity corresponding to 
the attractor configuration, one acting as an "atractor" 
to the dynamics and the other as a "repeller". With 
respect to the latter point, this system has mathemati- 
cally the characteristics of a chaotic scattering system 
with two abosolute outcomes consisting of (i) escape 
to infinity or (ii) recollapse to the singularity. The pres- 
ence of this critical point is responsible for a rich and 
complex dynamics, engendering in phase space topo- 
logical structures such as homoclinic orbits to a center 
manifold. The physical singularities are the main point 
in the whole discussion, that is, when any one of the 
scale factors crosses zero, meaning a recollapse of the 
universe. As showed in 111311 any homoclinic crossings 
present in the dynamics of the classical model is not 
seen by the physical world since the mandatory recur- 
rence is lost because the physical dynamics has to be 
restricted to A(t) and B(t) greater than zero. There- 
fore the only separation between recollapse and escape 
to the attractor at infinity are the unstable and stable 
manifolds corresponding to the center manifold asso- 
ciated to the Einstein singularity. This establishes the 
difference between physical and mathematical integra- 
bility: in spite of the chaotic dynamics present in the 
equations the physical meaningful region does not see 
it (see also Qj). 

In the present work we study the dynamics of the 



Bianchi IX cosmological model as in IU3I1 to which 
we add a term of quantum potential inspired by the 
work of Alvarenga et al. 11511 whose presence repre- 
sents exactly the short-range effects due to the quan- 
tum behavior of matter in small scales and plays the 
role of a repulsive force near the singularity. In this 
work a similar term has been introduced in an "ad 
hoc" manner. As it will be seen, this potential re- 
stricts the dynamics of the model to the positive val- 
ues of A(i) and B(t) and alters some qualitative and 
quantitative characteristics of the dynamics of the clas- 
sical model. We show the common features of a large 
class of such potentials which depends only on a so- 
called r-equivalent variable: (AB 2 ) 1 / 3 . Picking a par- 
ticular example, we make a complete analysis of the 
phase space of the model finding critical points, peri- 
odic orbits, stable/unstable manifolds using numerical 
techniques such as Poincare section, numerical contin- 
uation of orbits and numerical globalization of invari- 
ant manifolds. We compare the classical and the quan- 
tum models and verify that the addition of this quantum 
term allows the existence of homoclinic crossings of 
the stable and unstable manifolds in the physical mean- 
inful region of the phase space (both A(i) and B{t) 
positive) thus allowing chaotic escape to inflation as 
well as chaotic bouncing near the singularity due to a 
new center-center equilibrium point. 



II. THE AXISYMMETRIC BIANCHI IX MODELS 
WITH A QUANTUM POTENTIAL 

In this section we present the system and examine 
the main characteristics of the phase space of the ax- 
isymmetric Bianchi IX with a quantum potencial to- 
gether with the classical model so that the distinctions 
between the two models become apparent. We study 
the main characteristics of the phase space which can 
be obtained from linear analysis, such as the equilib- 
rium points and their nature; we present also the in- 
variant plane which is similar to the classical Bianchi 
IX model. 

Axisymmetric Bianchi-IX cosmological models are 
characterized by two scale functions A(t) and B{t) 
with the line element |l6tl 



ds 2 



dt 2 ~A 2 (t)(u 1 ) 2 -B 2 (t) x 
(uj 2 ) 2 + (oj 3 ) 2 ^ 



(1) 



Here t is the cosmological time and the to 1 are Bianchi 
IX 1-forms satisfying dof = e* J u^ A ui . The mat- 
ter content of the models is assumed to be a pressure- 
less perfect fluid, namely dust, with energy density 
p and four velocity field Sq in the comoving coordi- 
nate system used, plus a positive cosmological constant 
A. The dynamics of the scale factors A(t) and B(t) 
is given by Einstein's equations, which are equivalent 



to Hamilton's equations generated by the Hamiltonian 
constraint 



H 



PAPB 



A P 2 A A 3 



AB 8B 2 
- 2AAB 2 - E a = 



2A 



IB 2 



(2) 



where pa and ps are the momenta canonic ally conju- 
gated to A and B, respectively, and Eq is a constant 
proportional to the total matter content of the model, 
arising from the first integral of the Bianchi identity: 
pAB 2 = E Q . 

The quantum term that we introduce in this paper 
is inspired by the results Alvarenga et al. lfl5ll for 
the Friedmann-Robertson- Walker with a single scale 
factor representing the radius of the universes. These 
authors study the behavior of the wave functions de- 
rived from the Wheeler-DeWitt equations. Accord- 
ing to the causal Broglie-Bohm interpretation of quan- 
tum mechanics these wave functions are associated to 
particles whose trajectories approach classical trajec- 
tories in the classical limit. However in the regions of 
quantum validity the behavior differs from the classi- 
cal one. They found that this difference is equivalent 
to the addition of an ad hoc potencial term to the clas- 
sical model, equivalent to a repulsive force. This term 
agrees with the behavior of the trajectories in the quan- 
tum limit. In the work of Alvarenga et al. the quantum 
potential term has the following expression: 






(3) 



where a is the radius of the universes. Since our model 
has two scale factors A(t) and B(t) with volume of 
space given by AB 2 we have chosen an average radius 
given by r = (AB 2 ) 1 / 3 and write the new Hamiltonian 
as: 



H 



PAPB 



A P \ 



2A 



A 3 
2B 2 



AB 8B 2 
- 2AAB 2 + aV(r) - E = 



(4) 



where, as before, pa and ps are the momenta canon- 
ically conjugated to A and B, respectively, and a is a 
constant. We note that Hamiltonian (@]i with a = 
reduces to the classical one ©, However it is clear 
that the new Hamiltonian is not compatible with Ein- 
stein's equations and the quantum term is completely 
"ad hoc". We now compare the linear behavior and 
some other features of both the classical and the modi- 
fied Hamiltonian. 

For cr = 0, the classical case, the corresponding 
Hamilton's equations have an equilibrium point whose 
coordinates are determined by a very easy equation and 



E : pa = Pb = 0, A = B = 



4A 



(5) 



with associated energy Eq = E cr u = \/jt- The 



eigenvalues are A c i i2 = ±V / AA C 3 ! 4 = ±2iy/2A which 
characterize the equilibrium point E as a saddle-center. 
For the value of A = i chosen in this work we have 

A c i, 2 = ±V2, A c 3 ,4 = ±iV2 at A = B = 1. 

When the quantum term is added, the equation for 
the equilibrium points has to be determined numeri- 
cally for each pair of parameters (A, <x). Neverthe- 
less, we can make a few general statements about such 
potentials, which can be written as functions of r = 
{AB 2 ) 1 / 3 . To begin with, the divergence of the poten- 
tial at A = B = is a desirable feature to avoid col- 
lapse: the FRW singularity probably disappears. Also, 
the quantum system keeps the invariant plane at A = B 
and pa = Pb/2. The Hamiltonian describing the dy- 
namics in the invariant plane has then the following ex- 
pression: 



-i n .2 o/l 

H(A, PA ) = 2g- + ^ - 2AA 3 + aV(A) (6) 

with (7 = for the classical case. The dynamics 
in these planes are simply the level sets of the above 
Hamiltonian function. 

The fixed points on this invariant plane can also 
be easily described, regardless of the particular form 
of V(r). The introduction of the quantum potential 
slightly shifts the classical fixed point, Aq, towards 
smaller values and a new one, A\, appears between 
and Aq. While the former is always a saddle-center, the 
latter can be either a center-center, if 



d 2 V 



dA 2 



> 



ZA X 



(7) 



or a saddle-center otherwise. It is clear that an equality 
in the above equation would correspond to a bifurca- 
tion point. We recall that Aq is the classical saddle- 
center fixed point, given by ((5]). 

Moreover, for both cases, a straightforward analy- 
sis of the infinity of phase space shows that it has two 
attractors in this region. One acts as an attractor (sta- 
ble configuration) and the other as a repeller (unsta- 
ble configuration). The scale factors, A and B, in the 
classical case (er = 0) approach the attractor as A = 
B ~ exp {y/k/it) and p B = 2p A ~ exp (2^/A/3t) 
(cf. Figure 1). 

In order to further analyze the problem, we have to 
choose a especific form for V{r) and to pick particular 
values for the free parameters a and A. We make a 
direct generalisation of © and write 



V[r = (AB A ) 



2\l/3l = 



1 



(8) 



Taking the pair A = 1/4 and a — 0.01, which will 
be used throughout this work, we also find the saddle- 



Singularity 




FIG. 1: Phase space portrait of the invariant plane A — B, 
Pa = Pb/2. a) The separatrices S are characterized by the 
energy Eo — E cr it = l/v4A. The point A = B = is a 
degenerate singularity; b) The invariant plane for a 7^ 0: note 
that the quantum term removes the singularity A = B = 
and the trajectories bounce back in the vicinity of this point. 



center equilibrium point E , slightly dislocated (as ex- 
pected): 

E : pa=Pb= 0, A = B = 0.99317156, (9) 

with associated energy E = 1.001006820, higher than 
in the classical case. The eigenvalues for this equilib- 
rium point are A«i, 2 = ±0.49479869 and A 9 3 .4 = 
±1.42393682i, so that E continues to be as saddle- 
center. However, the inclusion of the quantum poten- 
tial generates another equilibrium point which is not a 
collapse given by 



E : p A =PB=0 7 A = B = 0.24194223, (10) 

and the corresponding eigenvalues are A 9 i2 = 
±2.4044353 and A« 3 ,4 = ±5.84525307i which char- 
acterize this equilibrium point as a center-center, in 
agreement with the criterium given at (0. This is the 
first significant change that the quantum term brings to 
the dynamics of the models. The corresponding conse- 
quences will be seen in Sec.lIVI 

As one can see in the invariant plane shown in 
Fig. (JTJ due to the new fixed point, the left unstable 



branch of the separatrix connects to the left stable one, 
creating a homoclinic loop and avoiding the collapse. 

In the vicinity of the saddle-center E it is always 
possible to find a set of canonical variables such that 
the Hamiltonian may be expressed as 11711 



remark that from ( fTTT > the center manifold of unstable 
peridic orbits is defined only for AE > 0. 



III. THE CENTER MANIFOLD OF UNSTABLE 
PERIODIC ORBITS 



H{qi,q 2 ,Pi,P2) 



O(3) + AE = 



2A(pi 



-«2) 

(11) 



where AE = (E C rit — Eq). The critical point E is 
located at the origin qi = q2 = Pi = P2 = 0, with 
Ecrit = Eq. Here 0(3) denotes higher-order terms in 
the expansion. In the approximation where C(3)-terms 
are neglected, the Hamiltonian is separable with the ap- 
proximate constant of motions given by the partial en- 
ergies E 2 = \/A(pl - g|)/2 and E 1 = V2A(pj + qf). 
Note that E\ is always positive. In a small neighbor- 
hood of E the motion will be the composition of rota- 
tional motion and hyperbolic motion connected respec- 
tively to the conserved quantities E\ and E 2 . To dis- 
play the topology of the linearized motion in a neigh- 
borhood of E we consider the two following possibili- 
ties. 

If E2 = and p^ = qi = the motion cor- 
responds to linear unstable periodic orbits te in the 
plane (qi,Pi). Such orbits depend continuously on the 
parameter Eq. The second possibility is E2 = and 
P2 = ±92, that defines the one-dimensional linear sta- 
ble Vs and linear unstable Vu manifolds. At the critical 
point E, these manifolds are tangent to the separatrices 
S of the invariant plane. The direct product of the peri- 
odic orbit te with Vs and Vu generates, in the linear 
neighborhood of E, the structure of stable (te X Vs) 
and unstable cylinders (te X Vu). Every orbit which 
constitutes these cylinders coalesce to the orbit te for 
times going to +00 or —00, respectively. The energy 
of any orbit on the cylinders is the same as that of the 
periodic orbit te ■ Two types of nonlinear extension 
can be considered for this separable linearized motion. 
The first is the nonlinear extension of the cylinders, 
away from the nonlinear periodic orbit, which will be 
the subject of Section IV. The second is the nonlinear 
extension of the rotational motion plane, that we start 
to discuss now. 

In fact, as shown in the next section, in the nonlinear 
regime (when higher order terms 0(3) of the Hamilto- 
nian are taken into account) the plane (qi,pi) of the ro- 
tational motion extends to a two-dimensional manifold, 
the center manifold, of unstable periodic orbits of the 
system 11811 . The intersection of the center manifold 
with the energy surface H(Eq) = is a periodic orbit 
from which a pair of cylinders emanates. The struc- 
ture of the center manifold of unstable periodic orbits 
and the associated cylinders for both Hamiltonians (O 
and (0| will be displayed in the next Sections, through 
a rather detailed numerical treatment. Finally we must 



In order to study the full dynamics it is necessary 
to first calculate the periodic orbits of the full sys- 
tem based on the preceeding linear analysis. These 
orbits can be continued to the periodic orbits of the 
full system by the numerical continuation technique of 
computing the family of te periodic orbits IU9I1 12011 . 
In general, to compute a T-periodic orbit of an au- 
tonomous system x = f{y), through xo, means to find 
the zeros of the function F(xq) = <p(T, xo) — xq = 0, 
where (p(t, xo) is the solution starting at xq when t = 
0. It is equivalent to dealing with the problem of find- 
ing a fixed point of a Poincare map 7 , (xo) of (p(t, xq) 
associated to a section £ transversal to the flow. Given 
£ and V(xq) as above we have now to find xq E £ such 
that 

T(x Q ) = (<p(t, x )nT,)-x = V{x )~x = 0. (12) 

Equation ( TT2T > can be transformed into a numerical con- 
tinuation problem which permits to follow the periodic 
orbits along a desired family. A convenient Poincare 
section is ps = 0, pb < 0. Throughout this paper 
we fix A = 1/4 and a = and a = 0.01, such that 
E crlt = 1.0 and E crit = 0.99317156 for the classi- 
cal quantum cases, respectively. We start to determine 
the ^o -parametric family (0 < Eq > 1) of non-linear 
periodic orbits continued from the linear approxima- 
tion around E. The predictor-corrector method uses 
the solution of the variational equations atthe trial so- 
lution taken from the linear approximation, reduced to 
the section ps = 0; and the differential of this Poincare 
mapping is checked for area preservation at each step. 

In Figure |2] (top) we show both the (A, pa) and 
(B,ps) proj ections for the quantum case for a = 0.01. 
The shape of the center manifold is similar in both 
cases (classical and quantum), and as in the classical 
case, the periodic orbits avoid B = 0. 

The above pictures show that the center manifold is 
a shell-like surface in phase space. From this shell the 
stable and unstable manifolds arrive/emanate asymp- 
totically. They act as guides for nearby trajectories 
either taking them to the neighborhood of the cen- 
ter manifold or away from this neighborhood towards 
other regions of the phase space. An important ques- 
tion we want to answer is whether these manifolds col- 
lapse and cross either A = or B = before being 
able to escape to the attractor like in the classical case 
discussed in 11131 11-411 . In order to see this we want to 
find where (in a convenient Poincare section) there is 
a first intersection of the stable and stable manifolds 
(if they indeed do cross) signalizing homoclinic points 



FIG. 2: Projection of the center manifold with a = 0.01: 
(top) (A, pa) plane; (bottom) (B,pb) plane. Note that the 
periodic orbits are qualitatively similar to classical ones. 



of this map which correspond to homoclinic orbits of 
the full dynamics. The existence of such points in the 
physically meaningful region will be discussed in the 
next section together with a summary of the methodol- 
ogy. In the classical case, such homoclinic points occur 
only in the non-physical region, where A < Ill3l[l4ll 
since B is never negative. 



IV. HOMOCLINIC INTERSECTIONS OF THE 
UNSTABLE AND STABLE MANIFOLDS 

The extension of the linear unstable/stable manifolds 
can cross each other creating homoclinic orbits and a 
chaotic region dynamics, and indeed this is the case 
as our subsequent numerical analysis shows. We will 
see that the addition of the quantum potential brings 
the stable and unstable manifolds to the the physically 
significant region as opposed to the classical case IU3I1 . 
This means that the dynamics of the homoclinic tangle 
can fully act on the question for the Bianchi IX models 
with an "ad hoc" quantum potential: how many times 
a trajectory bounces before escaping to the de Sitter 
region. 

We first (numerically) calculate the stable and unsta- 
ble cylinders by a technique denoted Numerical Glob- 
alization of Manifolds lfl9l 12011 . Second, in a Poincare 
section transverse to a periodic orbit (here, ps — 
0,Pb < ), so that the orbit is seen as a point and 




FIG. 3: (A, pa) projections of unstable cylinder manifolds: 
branch of the unstable manifold which goes towards A — : 
for (top) a = and (bottom) a — 0.01. Note that the ad- 
dition of the quantum potential makes the unstable manifold 
avoid the collapse by bouncing back towards the singularity. 



the cylinders are curves tangent to the directions of the 
eigenvectors of the Monodromy Matrix corresponding 
to the real eigenvalues. We note that these curves di- 
vide the energy space in interior and exterior regions 
and are formed by the iteractions of the orbits belong- 
ing to an invariant manifold. Therefore, orbits in each 
one of these two regions are subject to different out- 
comes under further iteractions if there are homoclinic 
intersections, as discussed later on. 

The symmetry of the problem in the section 
((A,pA,t) —> (A, —pa, —t)) allows saving computa- 
tional effort: from the Poincare map of the unstable 
manifold of a periodic orbit of the center manifold, we 
can obtain the corresponding Poincare map of the sta- 
ble one by just applying this symmetry. The branches 
of the stable and unstable manifolds which go towards 
the origin are the ones that intersect each other. The 
opposite branches go to the attractor-repeller/5o/«fi. 

In Figure [3] (bottom) we show the branch of the 
unstable manifold, in the quantum case (a = 0.01), 
for Eq = 0.97112431 which goes towards the region 
A — B — without collapsing, and returns again 
to a neighborhood of the periodic orbit from which it 
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FIG. 4: From top to bottom, (top) a = (classical case) four 
cuts (A, pa) of the unstable and unstable cylinder in the sur- 
face of section (p B = 0,p B < 0), for E = 0.9706090, 
0.90620828, .70789778 and 0.56639600; (bottom) four 
cuts as above for a = 0.01 for the E = 1.0095581, 
1.00706088, 1.00274454 and 0.99685469. Note that the 
quantum term brings the Poincare section of the tubes to the 
A > region while they are in the A < region in the 
classical case (see top figure). 



FIG. 6: (A, pa) (top) second and (bottom) third iteractions of 
the Poincare map in the surface of section (j>b = 0,Pb < 0) 
of orbits on the unstable manifolds to the p.o. with Eo — 
1.00706088. Note that the external shape of the second cut 
is practically the same as the first cut shown in Fig(5] The 
external part of the third cut mimics the second one. 




FIG. 5: (A, pa) first iteraction of the Poincare map in the 
surface of section (jps = 0,Pb < 0) of orbits on the 
unstable and stable cylinder manifolds to the p.o. with 
Eq — 1.00706088. Note that they are mirror images of 
each other with respect to pa = due to the symmetry 
(A, pa, t) •<— (A, —pa, —t) in the Poincare surface of sec- 
tion. The two intersections represent the 1-turn homoclinic 
orbits. 



FIG. 7: (A, pa) superposition of the first and second iterac- 
tions of the Poincare map in the surface of section (ps = 
0, pB < 0) of orbits on the unstable and stable manifolds to 
the p.o. with£ = 1.00706088. 




FIG. 8: The unstable manifold branch (top) at the onset of 
very large values of pa at Eq = 1.00274454; (bottom) for 
-Bo = 0.99585469, that is, for a larger value of energy the 
density of trajectories that approach A — increases, but 
part of them approach the periodic orbit and escape to the de 
Sitter region. 



emanated. Some orbits escape to the attractor-repeller 
points while the remaining ones return again towards 
A = B = 0. This is an indicative of homoclinic 
chaos. In Figure [3] (top) we show the corresponding 
branch of the unstable manifold for the classical case 
(<7 = 0). It emerges from a periodic orbit with en- 
ergy Eo = 0.99470278 and there are also escaping and 
returning trajectories. However, in the classical case 
this branch of the unstable manifold crosses the sin- 
gularity A — 0, indicating that the physical collapse 
takes place before it has time to return to the neigh- 
borhood of the periodic orbit. The same is true for the 
stabe manifold, which in this projection is mirror sym- 
metric with respect to the A — axis (not shown here 
to avoid cluttering the figures). So far the non collaps- 
ing branches of these manifolds are the main difference 
between classical and quantum models. Note that in 
the latter case, the singularity is avoided. This feature 
could have been guessed from the difference between 
the invariant planes shown in Figures Q] (top) and (bot- 
tom). 

In Figure [4] (top) and (bottom) we show four cuts 
of unstable manifolds in the surface of seccion (ps — 
0,Pb < 0) for the classical case (a = 0) and quan- 



tum case (a = 0.1), respectively, for decreasing val- 
ues of the energy as shown in the label of the Figure. 
The introduction of the quantum potential changes rad- 
ically the scenario: the cuts are brought to the physical 
significant region A > 0. We note that as the trajec- 
tories of the tubes going towards A — the momen- 
tum pa — > oo. However for energies near the criti- 
cal point the stretching effect is not so strong. In the 
classical case it was shown in 111311 that there are two 
homoclinic orbits that are seen in the Poincare section 
as the intersection of the two curves representing the 
unstable manifold which occurs at A < 0. As can 
be seen in Figure [5] where we show both the unstable 
and stable manifolds sections, there are two points of 
intersection: one near A = 0.14 and the other near 
A = 0.8. The energy of the periodic orbit, as well as 
of the stable/unstable manifolds is E = 1.00706088. 
Note that the sections of these two manifolds are mir- 
rors of each other as commented above. The two in- 
tersections represent homoclinic orbits comming from 
the periodic orbit through the unstable manifold, mak- 
ing one turn and returning assymtoptically to the peri- 
odic orbit through the stable manifold. Figures [6] (top) 
and (bottom) exhibt the second and third cuts of the 
unstable manifold, respectively. As can be seen in the 
second cut the set of orbits which do not escape form 
a very convoluted figure which gets larger and thinner 
as it tends to the first cut from its inside. The third cut 
which can be seen in Figure [6] (bottom) is even more 
convoluted since it has to goes towards the second cut: 
note that in this figure the "contour" of the first cut is 
denser due to the accumulation of points from the sec- 
ond and third cut. We also observe that the third cut is 
folded back inside the second one. Figure|7]shows two 
cuts of the unstable and stable manifolds and it is clear 
that higher-order homoclinic points have been created 
besides the first two ones shown in Figure[5] These or- 
bits are asymptotic to the periodic one after two turns. 
In this way a whole mesh of homoclinic points and 
initial conditions on different sides of these points can 
lead to escape after a given number of bounces. 

Therefore it is clear that the addition of the quan- 
tum potential creates a homoclinic tangle contained 
in the yl-and-B-greater-than-zero region, meaning that 
the escaping ruled by this homoclinic tangle has no 
chance to collapse. However, as the energy decreases 
the cuts go near A = and become more elongated 
(see Figure [8} and the pa moment becomes very large. 
This feature needs to be further examined. Note that 
the number of orbits having very large moment before 
turning towards the periodic orbit increases with de- 
creasing total energy. 



V. FINAL COMMENTS AND PERSPECTIVES 

We study one case of a class of quantum-inspired 
potentials in the classical axisymmetrical Bianchi IX 



models, namely V[r = (AB 2 ) 1 ' 3 ]. After pointing out 
the general features of such potentials, we focus on a 
particular expression, similar to the one exactly calcu- 
lated for FRW model Hal which is equivalent to adding 
a repulsive potential near the origin. The effect of this 
quantum term is to avoid the singularity and to give rise 
to an inflationary escape to the de Sitter region ruled by 
a chaotic scattering in opposition to the classical case 
where the two possible outcomes are separated only by 
the stable/unstable manifolds without the richness of 
the homoclinic tangle of these manifolds. 
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